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Abstract 

Consider the problem of scattering of electromagnetic waves by a doubly periodic struc- 
ture. The medium above the structure is assumed to be inhomogeneous characterized com- 
pletely by an index of refraction. Below the structure is a perfect conductor or an imperfect 
conductor partially coated with a dielectric. Having established the well-posedness of the 
direct problem by the variational approach, we prove the uniqueness of the inverse problem, 
that is, the unique determination of the doubly periodic grating with its physical property 
and the index of refraction from a knowledge of the scattered near field by a countably infi- 
nite number of incident quasi-periodic electromagnetic waves. A key ingredient in our proofs 
is a novel mixed reciprocity relation derived in this paper. 

Keywords: Uniqueness, Maxwell's equations, inhomogeneous medium, doubly periodic 
structure, mixed boundary conditions, mixed reciprocity relation, inverse problem. 

1 Introduction 

Scattering theory in periodic structures has many applications in micro-optics, radar imaging and 
non-destructive testing. We refer to |20| for historical remarks and details of these applications. 
This paper is concerned with direct and inverse problems of electromagnetic scattering by a 
doubly periodic structure. The medium above the structure is assumed to be inhomogeneous. 
Below the structure is a perfect conductor which may be partially coated with a dielectric. 
Let the doubly periodic structure be described by the doubly periodic surface 

Ta := {xGM 3 |x3 = /(xi,x 2 )}, 

where / G C 2 (IR 2 ) is a 27r-periodic function of x\ and x 2 - 

f(xi + 2ni7r, x 2 + 2n 2 ir) = f(x\, x 2 ) Vn = (ni,n 2 ) G I?. 

Assume that the medium above the structure T\ is filled with an inhomogeneous, isotropic, 
conducting or dielectric medium of electric permittivity e > 0, magnetic permeability fj, > 
and electric conductivity a > 0. Suppose the medium is non-magnetic, that is, the magnetic 



1 



permeability jjl is a fixed constant in the region above T\ and the field is source free. Then the 
electromagnetic wave propagation is governed by the time-harmonic Maxwell equations (with 
the time variation of the form e~ wt , oj > 0) 

curlE 1 - iujfiH = 0, curli? + uo(e + ia/uj)E = 0, 

where E and H are the electric and magnetic fields, respectively. Suppose the inhomogeneous 
medium is 27T-periodic with respect to the x\ and x 2 directions, that is, for all n = (m, n 2 ) G Z 2 , 

e(xi + 27rni,X2 + 27rn2,x 3 ) = e(xi,x 2 ,x 3 ), <j(xi + 2irni,x 2 + 27rn 2 , 0:3) = a(xi,x 2 ,x 3 ). 

Suppose above the structure Ti is another doubly periodic surface defined by 

T := {x G R 3 I x 3 = g(x 1 ,x 2 )}, 

where g G C 2 (R 2 ) is a 27r-periodic function of x\ and X2: 

5(0:1 + 2ni7r, a; 2 + 2n 2 7r) = g(x 1: x 2 ) Vn = (ni, n 2 ) G Z 2 , 

which separates the region above T\ into two parts: 

tt := {x G M 3 I x 3 > 5(0:1,0:2)}, 

fii := {x G M 3 I f(xi,x 2 ) < x 3 < g(xi,x 2 )}. 

Assume further that e(x) = €q, a = for x G f2o (which means that the medium above the layer 
is lossless) and that the doubly periodic surface Ti is a perfectly conductor coated partially with 
a dielectric. 

Consider the scattering of the electromagnetic plane wave 

E\x) = pe ikox - d , H\x) = re ikox - d 

incident on the doubly periodic structure Tq from the top region Qq, where ko = ^/eoftuj is the 
wave number, d = (ai,a 2 ,—f3) = (cos 9\ cos 6 2 , cos 9\ sin #2, — sm ^1) is the incident wave vector 
whose direction is specified by 9\ and 9 2 with < 9\ < it, < 9 2 < 2ir and the vectors p 
and r are polarization directions satisfying that p = 1/ ]I/eo(r x d) and rJ-d. The problem of 
scattering of time-harmonic electromagnetic waves in this model leads to the following problem 
(the magnetic field H is eliminated): 



curl curl E - k%E = 





in 


U , 


(1.1) 


curl curlS — klq(x)E = 





in 


fil, 


(1.2) 


v x curlE 1 1+ = Ao^ x curl_E 




on 


To, 


(1.3) 


v x E = 





on 


Tl,D, 


(1.4) 


v x cm\E — ip(v x E) x v = 





on 




(1.5) 



where q(x) = (e(x) + ia(x)/uj)/eo is the refractive index, v is the unit normal at the boundary, 
E = E % + E s is the total field in with E s being the scattered electric field, Ti = U Ti^, 
Ao and p are two positive constants. 

We require the scattered field E to be a -quasi-periodic with respect to x\ and x 2 in the sense 
that E(x\, x 2 , x 3 )e~ ia ' x is 2tt periodic with respect to x\ and x 2 , where a = (ai,a2,0) G M 3 . It 
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is further required that the scattered field E satisfies the following Rayleigh expansion radiation 
condition as X3 — > +00: 

E s {x) = V En e< a -- x+ ^\ x 3 > g+ := max 5 (x 1 , x 2 ), (1.6) 

' 2:1,2:2 

nez 2 

where ce n = («i + rai, 02 + n 2> 0) £l 3 , S n = (E^ ,En ,En) £ C 3 are the Rayleigh coefficients 
and 

= f (A; 2 - |a„| 2 ) 1/2 ifK| 2 <^ 

^ \z(K| 2 - ^) 1/2 if|«n| 2 >fc 2 

with z 2 = —1. From the fact that div E s (x) = it is clear that 

a n ■ E n + p n E^ = 0. 

Throughout this paper we assume that /3 n 7^ for all n £ Z 2 . 

The direct problem is to compute the scattered field E s in Qq and E in Q\ given the incident 
wave E l , the diffraction grating profiles Tq and T\ with the corresponding boundary conditions 
and the refractive index q(x). Our inverse problem is to determine the grating profile Ti together 
with the impedance coefficient p in the case when the interface grating profile Tq is known and 
the refractive index q in the case when the grating surfaces To and T\ are known and flat, utilizing 
the knowledge of the incident wave E 1 and the total tangential electric field v x E on a plane 
F/i = {x G K 3 I X3 = h} above the inhomogeneous layer. 

Problems of scattering of electromagnetic waves by a doubly periodic structure have been 
studied by many authors using both integral and variational methods. The reader is referred to, 
e.g. [U El HI [5J EH E21 [IS1 H] for results on existence, uniqueness, and numerical approximations 
of solutions to the direct problems. Compared with the direct problem, not much attention 
has been paid to inverse problems from doubly periodic structures although they are not only 
mathematically interesting but have many important applications. For the case when ri j = 
and the medium above the periodic structure T± = Ti t jj is homogeneous, the inverse scattering 
problem has been considered in [21 [71 E]. If the medium is lossy above the perfectly reflecting 
periodic structure, Ammari [2] proved a global uniqueness result for the inverse problem with one 
incident plane wave. If the medium is lossless above the perfectly reflecting periodic structure, 
a local uniqueness result was obtained in [7] for the inverse problem with one incident plane 
wave by establishing a lower bound of the first eigenvalue of the curl curl operator with the 
boundary condition (jl.4p in a bounded, smooth convex domain in 1R 3 . The stability of the inverse 
problem was also studied in [7J. Recently in [6], for the class of perfectly reflecting doubly 
periodic polyhedral structures global uniqueness results have been established in [6] for the 
inverse problem in the case of lossless medium above the structure, using only a minimal number 
(though unknown) of incident plane waves . Further, for a general Lipschitz, bi-periodic, partly 
coated structure Ti a global uniqueness result was proved in |13j for the inverse problem in the 
case of a lossless, homogeneous medium above the structure, using infinitely many incident dipole 
sources. 

On the other hand, for the case when T\j = (i.e. Ti = Tid), Ao = 1 and the grating 
surfaces Tq and T± are known and flat, a global uniqueness result was obtained in |14| for 
reconstructing the refractive index q, using all electric dipole incident waves (see [15] for the 
corresponding result in the 2D case). 
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In this paper, we prove global uniqueness results for the inverse problem of recovering a gen- 
eral smooth bi-periodic profile with a mixed boundary condition and a known bi-periodic interface 
from a knowledge of near field measurements above the known interface with a countably infi- 
nite number of quasi-periodic incident waves E l (x; m) = (1/fcg) curl curl [pexp(ia m ■ x — i^ m xs)], 
m = (mi, m-i) G 1?. Further, we also establish a global uniqueness result for the inverse problem 
of determining the refractive index q which depends on only one direction {x\ or X2) for the case 
when ri / = (i.e. T\ = T\ t o) and the grating surfaces Tq and Ti are known and flat, using a 
countably infinite number of quasi-periodic incident waves E l (x;m). This is an improvement to 
the result of |14j . A key ingredient in our proofs is a novel mixed reciprocity relation derived in 
this paper for bi-periodic structures. 

The rest of the paper is organized as follows. In Section [21 we introduce some suitable quasi- 
periodic function spaces and the Dirichlet-to-Neumann map on an artificial boundary above the 
structure. The problem (|l.ip - (|1.6p is then reduced to a boundary value problem in a truncated 
domain. In Section [3l we establish the well-posedness of the scattering problem (|l.ip - (jl.6p . 
employing a variational approach. Section |4] is devoted to the inverse problems. In Subsection 
14.11 novel mixed reciprocity relations are established for doubly periodic structures, which play 
a key role in the proofs of the uniqueness results for our inverse problems. Subsection 14.21 is 
devoted to the unique determination of the doubly periodic grating profile I\ with its physical 
property, where it is assumed that the interface Tq is known and the refractive index q is a known 
constant. Subsection 14.31 is concerned with the unique reconstruction of the refractive index g, 
where we only consider the case when the shape of the two grating profiles is known and flat, 
which improves the result in |14| . 



2 Quasi-periodic function spaces 

In this section we introduce some function spaces needed in the study of our problems. Due 
to the periodicity of the problem, the original problem can be reduced to a problem in a single 
periodic cell of the grating profiles. To this end and for the subsequent analysis, we use Tj,Qj 
(j = 0, 1), Ti^D,Tij and for h £M. again to denote the single periodic part (i.e. in the range 
< xi,x% < 2tt) of the corresponding notations defined in the last section. We also need the 
notation 

tt h = {x G M 3 I < xi, x 2 < 2ir, f(xi,x 2 ) < x 3 < h} 

for h > max{/(x') | x' G R 2 }. 

We now introduce some vector quasi-periodic Sobolev spaces. Let 

#(curl , U h ) = {E{x) = Y, E n (x 3 ) exp(ia n -x)\E€ (L 2 (n h )f, curl E G (L 2 (^)) 3 } 

nGZ 2 



with the norm 

2 



E \\H(cwl,(l h ) ~ W E \\L2(Q h ) + H curl ^llL2 



(fifc) 



Note that the a-quasi-periodic space -ff (curl , f2&) is a subset of the classical vector space 
H(curl , Qb) defined by 

M(curl , fi 6 ) = {E G (L 2 (n b ) f | curl£ G (L 2 {Q b )) 3 } 
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with the norm ||-E||ji( cur j ^ = \\E\ | 2 2 + Hcurl-EH 2 -^^)- Further, it was shown in [5J that 
if (curl , £lb) can be characterized as 

ff(curl,n 6 ) = {E G M(curl,fi 6 ) \e 2niai E(0,x 2 ,x 3 ) x e x = E(2ir,x 2 ,x 3 ) x a, 
e 2nia2 E(x u 0,x 3 ) x e 2 = E( Xl ,2n,x 3 ) x e 2 }, 

where ej = (1,0, 0) and e 2 = (0, 1, 0). 

To deal with the mixed boundary conditions (II, 4ft and (11. 5p . we introduce the subspace of 
ff(curl,n h ): 

X:={Ee H{cval,n h ) \ v x E\ V%G = 0, i/ x £|r M G L?(ri,j)} 
with the norm 



\E 



lx = ll^lllr(curi,n h ) + ll vx ^lli?crix) 



where L^Ti,/) = {E G (L 2 (ri j7 )) 3 | • E 1 = on r 1( j}. 
For x' = (xi, x 2 , /i) G r^, s£R define 

H t s (T h ) = {E(x') = E ^ exp(m n • x') \ E n G C 3 , e 3 ■ E = 0, 

n^iiWi0 = £(i + w a mi a <+°°} 

nez 2 

H|(div , r h ) = {£(x') = E ^ exp(ia n • x') | -E n G C 3 , e 3 ■ E = 0, 

nez 2 

ll^llfr.(div ) r h ) = E C 1 + l«-| 2 ) S d^| 2 + \ E n ■ «n! 2 ) < +00} 
nez 2 

H*(cw\ , r A ) = {E{x') = E exp(ia n • x') | £7 n G C 3 , e 3 ■ E = 0, 

nGZ 2 

ll^ll^Ccurl,^) = E t 1 + l«n| 2 ) S (l^n| 2 + |£„ X a n | 2 ) < +Oo} 
nGZ 2 

and write L|(r^) = H^iT^). We have the duality result: 

(fl/fdiv,^))^^— ^curl,^). 
Recalling the trace theorem on EI (curl , H^), we have 

ff t ~ 1/2 (div , F h ) = {e 3 x S|r h I E G (curl , Sl h )} 

— 1/2 

and that the trace mapping from If (curl , U^) to -ff t (div , T/J is continuous and surjective (see 
[8] and the references there). We also need the trace space Y(Tq) and its duality space Y(Tq)': 

Y(T ) = {/G^^lVro-ZG^/^ro)} 

Y(r y = {/GFr 1/2 (r )iv ro x f e H-^im, 
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where Vr denotes the surface gradient on Tq. Note that the trace space Y(Vo) can also be 
defined as follows (see |10) and |17[ p. 58-59]): 

y(T ) = {/ G (#~ 1/2 (r )) 3 | there exists E G #(curl , fi^Hi) 
with u x E = f on Fq, v x £ = on T/j} 

for h > g + , with norm 

||/|| y(ro) = inf{[|^|| ff(curIinhX iSx) I ^ e ffCcurl.JlhXHO, !/ x £| ro = /, v x ^| r ^ = 0}. 
We assume throughout this paper that q satisfies the following conditions: 
(Al) q G C^Hi) and = 1 for x G Sl ; 

(A2) Im [q(x)] > for all x G H x if T^/ / and Im [<?(x )] > for some x G Hi if i\/ = 0; 
(A3) Re [q(x)] > 7 for all x G Oi for some positive constant 7. 

3 The direct scattering problem 

In this section we will establish the solvability of the scattering problem (|l.ip - (|1.6|) . employing the 
variational method. To this end, we propose a variational formulation of the scattering problem 

in a truncated domain by introducing a transparent boundary condition on for h > g+. 

-~ 1/2 

Let x' = {x\,X2, h) G for h > g + . For E G H t (div , T^) with 



nSZ 2 

define the Dirichlt-to-Neumann map 1Z : H t 1,/2 (div , T/J — > H t 2 (curl , r\) by 

(KE)(x') = (e 3 x curl£) x e 3 on 

where E satisfying the Rayleigh expansion condition (|1.6p is the unique quasi-periodic solution 
to the problem 



The map 1Z is well-defined and can be used to replace the radiation condition (jl,6p on IV 
The scattering problem (ll.ip - (ll.6p can then be transformed into the following boundary value 
problem in the truncated domain Q^: 



where f% = — v x E l \r G Y (Tq), f 2 = — v x curli? l |r G Y(Tq)', / 3 = = and, for any vector 
i 7 , (F)t = iy x F) x v denotes its tangential component on a surface. 




curl curl E — k 2 E = for x% > h, 



uxE = E(x') onY h . 



curl curl E — k^E = in J)q, 

curl curl — /cQgE 1 = in Oi, 

1/ X .El-f — i/ X E\— = /1, z^x curlE 1 ^ — \qu x curli?|_ = f% on To, 

v x E = /3 on Ti^, ^ x curli? — i/>(z^ x E) x v = f± on r^j 

(curlE) T - 7e(e 3 x£)=0 on r\, 



(3.1) 
(3.2) 
(3.3) 
(3.4) 
(3.5) 
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Remark 3.1. In the case when k^q(x) = k\ is a constant and the incident field is given by 
the electric dipole source E % [x) = Gi(x,yo)r for yo £ an d r G R 3 (e.g. in the problem 
(|£2>(|£g]l of Lemma ED, we have fx = v x W\ v G F(r ), f 2 = A ^ x curl£*|r G Y(T )', 
h = -vx E^r^ G Y(F hD ), f A = -isx curl^| rijI + ip(y x F?) x u\ TlI G L 2 t (T 1}I ) in the 
problem (|3.ip - (|3.5p . where Y(Ti jj) is defined in the same way as Y{Tq) with Tq replaced by 
ri, D (see [13]). 



Define 



Y := {E€fl"(curl,J2i)nfr(curl,fi h \ni)|i/x£?|r lil3 =/ 3 , 

v x £| riI G L?(ri,j), v x £| + - 1/ x E\- = h on T } 



with the norm 



l^llr - ll s llH(curi,ni) + ll- £7 llH( CU ri I n h \ni) + W v x ^HiKr^)- 



Then the variational formulation for the problem (|3.ip - (|3.5p is given as follows: find E GY such 
that 

A(E, F) = B(F) VFeX, (3.6) 

where 

A(E,F) := A / (curl E • curl F - k%qE ■ F) dx 

+ / (curl£ • curl F - k$E ■ F)dx 
Jn h \Qi 

-i\ p / E T ■ F T ds - / TZ{u x E) ■ (u x F)ds, 
B{F) := I f 2 .F T ds + X [ h-F T ds. 

If fi G Y(To), then there exits Eq G ff (curl , fi/j\Oi) such that v x -Eb|r = /i> ^ x ^o|r h = 0. 
Similarly, for / 3 G Yfi^u) there exists a function E\ G if (curl , fi^) such that z/ x E\\y x d = fs, 
v x E\\? h = and ^ x £^i|ri j £ Lf (T\ j). Let E = E — Eq — E\, where 2?o is a function in 
satisfying that Eq\q } ^ = £o and -Eblci = 0. Then E £ X and the variational problem (|3.6D is 
equivalent to the problem: find E G X such that 

A(E, F) = B{F) VFeX, (3.7) 

where B{F) = B(F) - A(E ,F) - A{Ex,F). 

Theorem 3.2. Assume that the conditions (A1)-(A3) are satisfied. Then the problem (|3.1|) — 
(|3.5p has a unique solution E G Y for any f\ G Y(Tq), fi G Y(TqY } G l^r^o) and f± G 
L|(ri/). Furthermore, we have 

\\E\\y < C(\\fx\\ Y (To) + ll/2||y(r )' + ll/3||y(r li£) ) + 1 lx,»(r XjJ -))> 
where C is a positive constant depending only on 0,^- 
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Proof. It is enough to prove that the problem (|3,7p has a unique solution E E X with the 
required estimate. 

We first prove the uniqueness of solutions. To this end, let fj = 0, j = 1, 2, 3, 4 and let F = E 
in ([57711 . Then E) = 0, that is, 

A f (|curl£| 2 - k 2 q\E\ 2 )dx + [ _ (|curl£| 2 - k 2 \E\ 2 )dx 

-i\ p / \E T \ 2 ds - / Tl(y x E) ■ (y x £)ds = 0. 

Taking the imaginary part of the above equation and noting that the imaginary part of the last 
integral in the above equation is non-negative (see [131 Equation (16)]), we deduce that 

ifeg / lm(q)\E\ 2 dx + p \E T \ 2 ds<0. (3.8) 

If Tij = 0, then by (|3.8|) and the condition (A2) we have E = in a small ball B{x$] 5) C fii. 
By [91 Theorem 6] we have E E (i? 1 ^)) 3 . Thus, by the unique continuation principle (see 
[191 Theorem 2.3]) we have E = in Qi. This, together with the transmission condition (|3.3p 
and Holmgren's uniqueness theorem, implies that E = in 0,h\^i- If Ti / 7^ 0, then fj3.8j) and 
the boundary condition f|3.4|) yield that f x -E|r 17 = and v x eurl-Elr^ = 0. By unique 
continuation principle again we have E = in Again, from the transmission condition (|3.3p 
and Holmgren's uniqueness theorem it follows that E = in f2/j\Ox- The uniqueness of solutions 
is thus proved for both cases. 

Now, arguing similarly as in the proof of Theorem 4.1 in |14j or Theorem 3.1 in |13j (see 
[14^ [T3] for details) we can prove that the problem (|3.7p has a solution E E X satisfying the 
estimate 

\\E\\ X < C(||/ 2 ||y(r )' + \\U\\ L 2(r ltI ) + H-^llirfouri.nAni) + IM*) 

with C depending only on Q^. Since E = E — Eq — E\, and by taking the infimum over all 
Eq E .ff (curl , r2/i\f2i) such that v x -Eb|r = /1 and v x i?o|r h = and over all E\ E .ff (curl , fi^) 
such that 1/ X £?i|r 1 ^ = /3, z/ X -Ei|r h =0 and v X E\\r 1 1 E L|(ri /) the desired estimate follows 
(on taking into account the definition of the norm on Y(Tq) and Y(Ti £>)). 

□ 



4 The inverse problems 

In this section we consider the inverse problems of determining the doubly periodic grating 
profile / with its physical property and the refractive index q from a knowledge of the incident 
and scattered fields. To this end, we need the free-space quasi-periodic Green's function 

Go(x,y) = — 5 y~* — exp(ia n • (x - y) +if3 n \x 3 - y 3 \) 

ngZ 2 Hn 

provided f3 n 7^ for all n E Z 2 (see |18j). In the neighborhood of x = y, Gq can be represented 
in the form Go(x,y) = $>(x,y) + a{x — y), where = exp(iko\x — y\)/(4w\x — y\) is the 
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fundamental solution to the three-dimensional Helmholtz equation (A + kfyu = and a(x — y) 
is a C°° function (see [16J for the 2D case). We now introduce the quasi-periodic Green's tensor 
Gq £ C 3x3 for the time- harmonic Maxwell equations: 



G (x,y) = G {x,y)I + -^V x dw x (G {x,y)I), x^y, (4.1) 

Kr. 



where I is a 3 x 3 identity matrix. Consider the following incident dipole source located at z £ R 3 
with polarization p (\p\ = 1): 

E l (x) = Gq(x, z)p, i/z. 

Clearly, we have 

curl curlS* (x) - k%E l (x) = 0, x ^ z. 
4.1 Mixed reciprocity relations 

We establish two mixed reciprocity relations for the doubly periodic structure, which play a key 
role in the proofs of uniqueness results for the inverse problems. 

Lemma 4.1. Assume that k$q(x) = kf is a constant. For m = (mi, 7712) £ let E l (x;m) = 
(l/fcg) curl curl \pex.p(ia m • x — i/3 m X3)] and let E(x\m) (which is the sum E l (x;m) + E s {x;m) 
in £1$) be the solution to the the scattering problem (|l.ip — (|1.6p with E' l (x) = E' l (x;m). On 
the other hand, define a := —a and for yo G £li and r 6 R 3 let E l (x;yo) = Gi(x, yo)r and let 
E(x; yo) solve the scattering problem: 

curl curl E — k^E = 
curl curlS — k\E = 
v x E\ + = v x S|_, v x curlS| + = Ao^ x curlSl- 
v x E = on v x curlE 1 — ip(v x E) x v = 

E(x;y ) = E^x-^q) + E s (x;y ) 
E(x; yo) = ^ E n (y ) exp(zS n • x + /3 

ngZ 2 

S n • E n + ^ n • E^ = 0. 

i/ere, S ra and f3 n are defined by 



in 


f2 , 


(4.2) 


in 


^i\{yo}, 


(4.3) 


on 


r , 


(4.4) 


on 




(4.5) 


in 


^i\{yo}, 


(4.6) 


in 


X3 > g+, 


(4.7) 






(4.8) 



S n = (-ai + ni, -a 2 + n 2 ,0) and /3 n 



and Gi(x,yo) * s defined by (|4.1|) iwi/i a n and A;q replaced by a n and k\, respectively. Then we 
have 

r ■ E(y ; m) = ^^/3_ m S_ m (y ) ■ P- (4.9) 
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Proof. Note first that E(x; m) and E(x; yo) are well-defined by the well-posedness of the di- 
rect scattering problem. Applying Green's theorem in Oi \ B(yo,5) and using the fact that 
contributions of the vertical line integrals cancel out due to the periodicity, we have 

= f {curlcurl £ (x; m) .[G 1 (x, ! ,„ W - E (x ; „).c„rlcu I l[e i (x, I/ „) r ]} <fa 

|z/ x cw\E(x) ■ [Gi(x,yo)r] — v x curl [Gi(x, yo)r] ■ E(x)\ ds 

(4.10) 



+ 

r JTi JdB s (y ) 

h+h + h 



where -B(yo, 5) is a small ball centered at yo with the radius 5 such that B(yo, 5) cOi. 

We now analyze the asymptotic behavior of I3 as 5 — > 0. From the definition that G\(x, yo) 
G\(x, yo)I + ki V x div x (G\(x, yo)I) it follows that 



h 



dB(y ,S) 



+ 

ldB(y ,S) 

h + h- 



v x curl E(x; m) ■ k 1 2 Vdiv [G\ (x, yo)r] — v x curl [G\ (x, yo)r] ■ E(x; m) 
v x curl£'(x; m) ■ rG\(x, yo)ds 



ds 



(4.11) 



The regularity of E(x; m) and the singularity of G±(x, yo) at x = Uo imply that I5 — > as 5 — > 0. 
On the other hand, by the divergence theorem on dB(yo, 5) it can be seen that 



/ 



dB s (y ) 
f 

9B s (y ) 
f 

9B s {yo) L 



v x curl E(x; m) ■ — ^-Grad div \G\ (x, yo)r] — v x curl \G\ (x, yo)r] • E(x; m) 

K 

Div(— v x curl E(x; m))— jdiv [Gi(x, yo)r] — v x curl [Gi(x, yo)r] ■ E(x; m) 

^1 



ds 
ds 



Idi 
Idi 



(y ■ curl curl E(x; m))-p>div [G\(x, yo)r] — v x curl [G\(x, yo)r] ■ E(x; m) 
kf 



ds 



(u ■ E(x; m))div [G±(x, yo)r] + v x E(x\ in) ■ curl [G±(x, yo)r] 



ds 



dB s (y ) 
r ■ E(y ;m) 



[v ■ E(x;m))VGi{x,y ) - VGi(x,y ) x (z^ x E(x;m)) 



ds ■ r 



as 5 — V 0. This combined with (|4.10|) and (14. lip implies that 
r ■ E(y ;m) 

v x curl E(x) ■ [G\(x, yo)r] — v x curl [G±(x, yo)r] ■ E(x) ds. (4.12) 



T JT U 

Similarly, we have on noting the regularity of E s (x; yo) that 



To .Ti 



v X curl E(x;m) ■ E s (x;yo) — u X curl E s (x;yo) ■ E(x;m) ds = (4.13) 
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Combine (jl7T2]) with (|3~T3j) to conclude that 
r • E(y ;m) 

v x curl E(x; m) ■ E(x; yo) — v X curl J5(x; yo) ■ to) 



r 



(is. 



Making use of the boundary conditions on Tj (j = 0, 1) and Green's theorem in f^f^i we obtain 
that 



r ■ E(y ;m) 

v x curl E(x; m)|_ • E(x; yo)\- — v x curl£'(x; yo) I— • E(x; m)\. 



r 
1 

-/ 

A Jr h 



ds 



v x curl E{x\ m)|+ ■ S(z; yo)|+ — ^ x curl£'(x; yo)| + ' to)|+ 



^ x curl m) • £/(ic; yo) — ^ x curl yo) • E{x; m) 



ds 



Now, by the Rayleigh expansion radiation condition and the divergence- free property for E s (x; m) 
and E(x;yo) we have, on noting that (3 n (a) = /3_ n (2), that 



v x curl£ ,s (x; m) ■ E(x; yo) — v x curlS(x; yo) " E s [x\ m) 



ds = 0. 



This implies that 



r ■ E(y ;m) 



A 7r h 

Insert E l (x;m) = A;Q~ 2 curlcurl [y ex.p(ia m 
ifinXz} into the above equation to get 



(is. 



^ x curl£'*(x; m) • yo) — ^ x curl yo) • E l (x; m) 

(3 m x 3 )\ and E(x;y ) = ^ £ n exp{iS n • x + 



r-E(y ;m) = y ^ {[^{.Vo) X e 3 ] x (a m ;-/3 m ) - e 3 x [(ct„;/3 n ) x E n (y )]^ e l 

/■27T /*27T 

'Pm / / e i(a " +Qm) :E dxidX2, 

Jo Jo 

where (a; 6) is defined as (a; 6) := a+ (0,0,6) and p m = p - [(a m ; — /3 m ) ■ p/ko\(a m ; — /3 m ). 

Finally, use the fact that a n + a m = (n + m, 0), /3-i(a) = /3_/(a) for all I £ Z 2 and 
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f*27T f'2lT 

/ e i{n+mfi) ' x dxidx 2 = for n + m / (0, 0) to conclude that 



'o J o 



r • E(y ;m) 



47T 2 i 



|[J5_ m (y ) x e 3 ] x (a m ;-/3 m ) - e 3 x [(a_ m ;/?-m) x £L m (y )]} • p m 

'Pm 



— <^ (a m ; -/3 m ) x [£_ m (y ) x e 3 ] + e 3 x [(a_ m ; /3_ m ) x E_ m (y )] \ 
o L J 

47r 2 i ( ^ ^ ^ ^ ^ 1 

|-/3 m £'_ m - [(a m ;-/3 m ) • .E_ m ]e3 + £l^(a_ m ; /3_ m ) -/3_ m £'_ m | -p r , 



A, 



|-Ei 3 i(yo)(2-m;/3-m) - 2/L m JL m (y )} -Pr, 



l<K 2 i 



Aq 

/3_m£-m(yo) -p, 



Ao 

where we have used the fact that 

(cK— mi fi—m) ' Pm — (o^mj /3m) ' Pm — 0) 
(^mi Pm) ' E— m — ( CX— m , /3_ m ) • E— m — 0. 

This completes the proof. □ 

If yo G Oo, define the total field E(x;yo) = E s (x;yo) + Go(x,yo)s in Qq, where Go(x,yo) is 
an a-quasi-periodic Green tensor defined in (|4.ip with a replaced with 3. Then arguing similarly 
as in the proof of Lemma l4.1l we can prove the following result. 

Lemma 4.2. Form = (mi, 7712) G Z 2 let E l (x;m) = (l/A^curlcurl [pexp(ia m • ^ — i/3 m x 3 )] and 
let E(x;m) (which is the sum E l {x;m) + E s (x;m) in $1q) be the solution to the the scattering 
problem (jl.ip — (|1.6p with E l (x) = E % (x;m). For y$ £ ^0, 3 = —a and r£l 3 Zei E l (x;yo) = 
Go(x,yo)r and let E(x;yo) (which equals to the sum E l (x;yo) + E s (x;yo) in ^io\{yo}) satisfy the 
Maxwell equations curl curl E — KqE = in &o\{yo} ara ^ curl curl E — k^qE = in fii together 
with the transmission condition 

v x + = ^ x f x curli£| + = Ao^ x curli?|_ on To, 

the boundary condition 

v x E = on r^o, ^ x cnvXE — ip(y x E) x u = on Tij 
and the Rayleigh expansion radiation condition 

E s {x;y ) = ^2 E n (y ) ex.p(ia n ■ x + /3 n x 3 ) for x 3 > g + , 

where 

a n -E n + X- E^ ] = 0. 

Then we have 

r ■ E{y ; m) = 8tt 2 i^ m E. m (y ) ■ p. (4.14) 
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4.2 Unique determination of the impenetrable profile / 

We now consider the unique determination of the impenetrable grating profile /, assuming that 
the interface profile g is known and k^qix) = k\ is a constant. A key ingredient in our proof is 
the mixed reciprocity relation for the doubly periodic structure (see Lemma l4.ip . 

Theorem 4.3. Assume that f3 n ^ for all n G Z 2 , the interface profile g is known and 
k$q(x) = k\ is a constant. Let /i, /2 £ C 2 (IR 2 ) be 2tt -periodic, let pi, p2 be two constants 
and let h > max a , eR 2{/i(x), /2(x)}. // v x Ef m \r h = v X -Ef mlr\ for all incident waves 
Em( x ) = (l/^o) cur ^ cur U e i exp(ia m • x — if5 m xs)] with m£Z 2 and I = 1,2,3, then 

fi = /2, T /lijD = Tf 2iD , T fuI = T h j, pi = p 2 , 

where ei is the unit vector in the direction xi, I = 1,2,3. Here, Ej >rn = E l m + Ej m in £Iq and 
Ej t m in ^lfj are the unique quasi-periodic solution of the scattering problem (jl.ip — (jl.6p with 
E % = E l m , p = pj and f = f j; where U fj = {x £ E 3 \ f j (x 1 ,x 2 ) < x 3 < g(xi,x 2 )}, j = 1, 2. 

Proof. We assume without loss of generality that f% ^ f2 and there exists a z* = (z\, z^,z^) £ Tfa 
with fi(z^,Z2) > f2(z\,z£), where Tf j = {x £ M 3 | x% = fj(xi,X2)}. We choose e > such that 
z e := z* + ee 3 £ Uf t fl Sl/ 2 . 

Let E e j be the unique quasi-periodic solution to the scattered problem (|4.2p - (|4.8p with yo = 
z e , p = pj, f = fj. By Lemma 14.11 we have 

r • E lm (z e ) = ^-^j3- m Ei,- m (ze) • e h (4.15) 
r ■ E 2i m(Ze) = -^-P- m E 2 - m {zt) ■ e h (4.16) 

where Ej tn (z t ) are the Rayleigh coefficients for E t j. 

On the other hand, from the Rayleigh expansion radiation condition and the assumption that 
v x Ef m \r h = v x E\ m |r ; j we conclude by the unique continuation principle that Ef m = i?! m in 
Qq- This, together with the transmission condition on Tq and Holmgren's uniqueness theorem, 
implies that E\ )in = i?2,m i n ^/i H Qf 2 i so E\ )m {z e ) = E2 t m(z t ). It then follows from (|4,15p and 
P~T6]) that 

— — P- m E\- m {z € ) ■ ei = — — P- m E2- m {z e ) ■ ei or ^ - m (.z e ) = E 2 ,- m (^e)- 

AO ^0 

Thus, by the Rayleigh expansion radiation condition we have E 6j i(x) = E e ^(x) for xz > g+- By 
the unique continuation principle, the transmission condition on To and Holmgren's uniqueness 
theorem again we obtain that 



E Et i(x) = E^ 2 {x) in £7 and E s e>1 {x) = E s e2 {x) in Q^nfl/ 2 . 

Without loss of generality we may assume that z* lies on the coated part of T^. Since z* has 
a positive distance from Tf 2 , then the well-posedness of the direct problem implies that there 
exists C > (independent of e) such that 

\(v x curli^ - ipiu x E s e l x v)(z*)\ = \(u x curl^ 2 - ip x v x E*^ x ^)(z*)| < C. 
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However, from the boundary condition on it is seen that 

\(u x curll^ - ipiv x E*^ x v){z*)\ 

= \(u x curl [Gi(-, z e )r] — ip\v x [Gi(-, z e )r] x v)(z*)\ — > +00 

as e — t- 0. This is a contradiction, which implies that f\ = /2, that is, VLf x = VLf 2 and = Tf 2 . 
Hence, we have Ei m = I?2,m m ^/i- We claim that r^^HT^/ must be empty (so Vt x D = ^f 2 ,D 
and r f 1: i = r f 2 j) since, otherwise, a similar argument as below deduces that the total field E\ jTn 
vanishes in fiy 15 which is impossible. 

Now let / = fx = /2- Then by the boundary condition we deduce that 

i(pi ~ Pi)(v x Ei >m ) x v = on 

If pi ^ P2, then the above equation implies that v x E\ m = on T\j, so by the boundary 
condition again v x curli?i jm = on Y\i. Thus, by Holmgren's uniqueness theorem, E\ m = 
in Oi. By the transmission condition on To and Holmgren's uniqueness theorem again it follows 
that -Ei,m = E l m + E( m = in Qq, which is a contradiction. The proof is thus completed. □ 

4.3 Unique determination of the refractive index 

We now consider the inverse problem of recovering the refractive index q. We only consider the 
case that T%j = 0, that is, the grating surface Ti is a perfect conductor. However, we expect the 
result to hold in a more general case by constructing special solutions of the Maxwell equations. 
Throughout this section we assume that the transmission constant Ao is known and the shape 
of the grating surfaces Tq and T± is also known and flat, that is, for two known constants b > c, 
g(x') = b and f{x') = c for all x' G M 2 . 

We have the following global uniqueness result for the inverse problem. 

Theorem 4.4. Assume that q = qj satisfies the conditions (Al) — (A3) and that qj depends 
only on x\ or x^, j = 1, 2. Let h > b. If 

v x E tm\r h = ^ x ElJ Th 

for all incident waves E l m (x) = (1//cq ) curl curl [e; exp(za m • x — if3 m xs)] with m S Z 2 and I = 
1,2,3, then we have q\ = qi- Here, Ej tTn = E % m + E s - m in VIq and Ej )Tn in £li are the unique 
quasi-periodic solution of the scattering problem (II. lh — (|1.6j) with E l = E % m and q = qj, j = 1, 2. 

Remark 4.5. Theorem 14.41 improves the result in [14, Theorem 5.4], where only the special case 
Ao = 1 is considered and incident waves of the form (I4.17P below are used for all r € L 2 (r^). 

To prove Theorem 14.41 we need the following denseness result which is related to the incident 
waves of the form 

E i (x;r)= [ G (x,y)r(y)ds(y), x 3 < h, (4.17) 

where r G L^(Th). This result was proved in |141 Lemma 5.2] for the case Ao = 1, and the general 
case can be proved similarly (see the proof of Lemma 5.2 in |14j). 



14 



Lemma 4.6. The operator F has a dense range in H t ^(div , To). Here, F : L^(Th) — > 

— 1/2 ^ ^ 

H t (div,ro) is defined by (Fr)(x) = e% x E(x;r)\- on Fq, where E(x;r) is the solution 
of the scattering problem (jl.ip — 1 1 .6|) with the incident wave E l (x) = E l (x;r) given by (|4.17p . 

Proof of Theorem 14.41 For any r £ (T^) and y G we have by Lemma FOl that 

r(y)-E°(y;m)=8ir 2 if3_ m E^ m (y)-e l , j = 1,2, I = 1,2,3, (4.18) 

where Ej_ rn {y) are the Rayleigh coefficients of the scattered field E s -{-;y) corresponding to q = qj 
and the incident wave E l (x) = Go{x,y)r(y). It follows from (|4,18p that 



r(y) ■ E s j (y;m)ds(y) = 8^/3_ m / _ m (i/)ds(i/) ■ e h (4.19) 

Jr h 

Denote by E?(x; r) and Ej(x; r) the scattered and total electric fields, respectively, corresponding 
to q = qj and the incident wave E l (x) = E l (x;r), j = 1,2. Then from the definition (|4.17p of 
E l (x;r) it is seen that 

/ Ej^ m (y)ds{y) are the Rayleigh coefficients of Ej(x;r). (4.20) 

On the other hand, from the Rayleigh expansion radiation condition and the assumption that 
v x Ef(x; m) = v x E^{x\ m) on we conclude on using the unique continuation principle that 
Ef(x;m) = E^x^m) in Q,q. This, together with (|4.19p and (|4,20p . implies that 



Ei.- m {y)ds{y) = I E 2 - m (y)ds(y). 

From this, the Rayleigh expansion radiation condition and the unique continuation principle it 
follows that 

Ef(x;r) = E 2 (x;r) or E\(x; r) = E 2 (x; r) in S7^\S7i. 

With the help of the transmission conditions on To, we get 

v X Ei(x;r)\^ = v x E 2 (x; r)|_ on r , 

v x curl E\(x\ r)|_ = v x c\xr\E 2 {x\ r)|_ on Tq. 

Now define E{x) := Ei(x;r) — E 2 {x;r) in Oi. Then E satisfies the equation 

curlcurlE 1 — k\q 2 E = k^(qi — q 2 )Ei(x; r) in Vt\ 

and the boundary conditions 

v x E = 0, v x curlE 1 = on To, 
i/xB = on Tj. 
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Thus it follows from Green's vector formula that 



<•■: 



,2 





(qi - q 2 )E 1 (x;r) ■ E 2 {x)dx 



f 
I 

Jn 



(curl curl E — /cqQ^I?) • E2(x)dx 



(curl curl E2 (x) — k^q^E^x)) ■ E(x)dx 







(4.21) 



for any r 6 Lf(Th), where E% G i/(curl,f2i) satisfies the Maxwell equation (jl.2p with q = q 2 
and the boundary condition v x -E^Ti = 0. 

Now by Lemma [4. 6 1 and ([4.21|) we obtain that 



where E\ satisfies of the Maxwell equation (jl,2p with q = qi and the boundary condition v x 



Finally, using the orthogonal relation f|4.22|) and arguing in exactly the same way as in the 
proof of Theorem 5.4 in |14| . we can easily prove that q\ = f/2- The proof is thus completed. □ 
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